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We study the improvement of simulations of QCD with dynamical Wilson fermions by combining the Hybrid 
Monte Carlo algorithm with parallel tempering. As an indicator for decorrelation we use the topological charge. 



1. INTRODUCTION 

Decorrelation of the topological charge in Hy- 
brid Monte Carlo (HMC) simulations of QCD 
with dynamical fermions is a long standing prob- 
lem. For staggered fermions an insufficient tun- 
neling rate of the topological charge Qt has been 
observed 610 . For Wilson fermions the tunneling 
rate is adequate in many cases |3|,|]]. However on 
large lattices and for large values of k near the chi- 
ral limit the distribution of Q t is not symmetric 
even after more than 3000 trajectories (see Fig- 
ure 1 of H and similar observations by CP-PACS 

The idea of parallel tempering is to improve 
transitions in parameter regions where tunneling 
is suppressed by opening ways through parame- 
ter regions with little suppression. In QCD the 
method has been applied successfully for stag- 
gered fermions ||. In (]]] parallel tempering has 
been used to simulate QCD with 0(a)-improved 
Wilson fermions without finding any gain, how- 
ever, with only two ensembles which does not take 
advantage of the main idea of the method. 

Here we use parallel tempering in conjunction 
with HMC to simulate QCD with (standard) Wil- 
son fermions in a parameter range close to the 
critical region, where the method is expected to 
be most advantageous. The gain achieved is 
demonstrated by studying time series and his- 
tograms of the topological charge. 



2. PARALLEL TEMPERING 

In standard Monte Carlo simulations one deals 
with one parameter set A and generates a se- 
quence of configurations C. The set A here in- 
cludes (3, k, the leapfrog time step and the num- 
ber of time steps. C comprises the gauge field 
and the pseudo fermion field. 

In the parallel tempering approach ||j9| one 
simulates N ensembles (Aj;Cj), i — 1, . . . , N in 
a single combined run. Two steps alternate: (1) 
update of N configurations in the standard way, 
(2) exchange of configurations by swapping pairs. 
Swapping of a pair of configurations means 



((AijCO.CAjjCj)) 



((Xi;Cj), (\j-,Ci)), if acc. 
((A*; Ci), (Xj-,Cj)), else 



with the Metropolis acceptance condition 



-Pawap(i) j) = min (l, 



(1) 



*Talk presented by H. Stiiben at Lattice '99, Pisa, Italy. 



AH = H Xi (Q) + H Xj (Cj) - H Xi (Cj) - H X] (Q). 

Since after swapping both ensembles remain in 
equilibrium, the swapping sequence can be freely 
chosen. Swapping neighboring pairs and proceed- 
ing in parameter direction towards criticality has 
turned out to be most advantageous. 

3. SIMULATION DETAILS 

We used the standard Wilson action for the 
gauge and the fermion fields and worked on an 8 4 
lattice. Our HMC program applied the standard 
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Figure 1. Comparison of standard with tempered HMC. The tempering data shown are from the run 
with 6 ensembles in Table 1. 



conjugate gradient inverter with even/odd pre- 
conditioning. The trajectory length was always 1. 
The time steps were adjusted to get acceptance 
rates of about 70%. In all cases we generated 1000 
trajectories (plus 50-100 trajectories for thermal- 
ization). 

Qt was measured by the field-theoretic method 
after 50 cooling steps of Cabibbo-Marinari type. 

4. RESULTS 

We have run several tempered HMC simula- 
tions in the quenched approximation (tempering 
in (3) and with dynamical fermions (tempering in 
k, at fixed (3 = 5.5 and f3 = 5.6). For compar- 
ison also standard HMC simulations have been 
performed. 

Figure 1 shows a typical comparison of results. 
We find that with tempering considerably more 
topologically nontrivial configurations occur and 
the histograms of Q t become more symmetrical. 



In standard runs Qt frequently stayed for quite 
some time near 1 or near —1, while with tem- 
pering this never occurred. The standard run 
at k — 0.156 shown in Figure 1, where Qt gets 
trapped in this way for about 200 trajectories, 
provides an example of this. Such observations 
have also been made on large lattices ||[| . 

While a correlation analysis cannot be carried 
out with the size of our samples, some quantita- 
tive account of the improvement by tempering is 
possible using the mean of the absolute change of 
Q t , called mobility in Q, 

^^Eiftw-^" 1 )!- ( 2 ) 

We give results for D\ in Table 1. The k- values 
of our tempered runs include the ones used by 
SESAM The considerable gain by tempering 
is obvious from Figure 1 and Table 1. 

For comparisons in Table 1 one should, in ad- 
dition to the individual errors at each of the 
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Table 1 

Mobilities D\ (sec (0)) on the 8 4 lattice at j3 = 5.6. The swap acceptance rates ([[]) achieved were about 
82% for Ak = 0.00025 and about 63% for Ak = 0.0005. 
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k- values (which are not small), take the over- 
all behavior within the particular runs into ac- 
count. Comparing the tempering results for 7 and 
6 ensembles it appears that adding intermediate 
points (not shown in Table 1) in this case leads 
to less improvement than extending the range to 
lower k- values does. The comparison of the re- 
sults for 21 and 6 ensembles possibly indicates 
that extending the range to very large K-values 
(with no transitions in standard simulations) af- 
fects the results at lower k. In any case our data 
show that by a clever choice of range and dis- 
tances of the parameter values the number of pa- 
rameter points can be strongly reduced. 

5. CONCLUSIONS 

Parallel tempering considerably enhances tun- 
neling between different sectors of topological 
charge and generates samples with more symmet- 
rical charge distributions than can be obtained by 
standard HMC. The histograms also get slightly 
broader or even become nontrivial thanks to this 
technique. 

The enhancement of tunneling indicates an im- 
provement of decorrelation also for other observ- 
ables. More satisfactory histograms are impor- 
tant for topologically sensitive quantities. Both 
of these features make parallel tempering an at- 
tractive method for large-scale QCD simulations. 
The method is particularly economical when sev- 



eral parameter values have to be studied anyway. 

A potential problem which remains to be stud- 
ied is that for a given parameter set the swap ac- 
ceptance rate (|l|) decreases for increasing lattice 
volume 0. In this context possibly the choice of 
range and distances of the parameter values gets 
more important. We therefore plan to extend our 
investigations to larger lattices. 
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